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Abstract. Due to computational requirements and numerical difficulties associated 

with coordinate singularity in spherical geometry, fully dynamic 3D magnetohydrody- 

namic (MHD) simulations of massive star winds are not readily available. Here we 

report preliminary results of the first such a 3D simulation using 6^ Ori C (05.5 V) 

as a model. The oblique magnetic rotator 6^ Ori C is a source of hard X-ray emitting 

Q/. , plasma in its circumstellar environment. Our numerical model can explain both the 

^^ ■ hardness and the location of the X-ray emission from this star confirming that magnet- 

• I ically confined wind shock (MCWS) is the dominating mechanism for hard Xrays in 

^5 ■ some massive stars. 

o 

,^' 1. Introduction 



Since the first detection by Einstein of X-ray emission from hot, luminous stars (iHamden et al. 
^ I 1 19791 ). a favored scenario has been that these X-rays arise from embedded shocks 

^ . gene rated by intrinsic instabilities in the radiatively driven stellar winds of such stars 



(Luc y & Whit3ll980l : lOwocki et alJ [T988: Feldmeier et al. 1997). For several O-type 



supergiants (e.g. Zeta Pup, Zeta Ori), the relatively soft, broadened emission line spec- 
(T^ \ trum observed by Chandra seem generally to support this paradigm. However, in sev- 

O ■ eral other hot stars (e.g. r Sco, 6^ Ori C, cr Ori E), the much harder, narrow emission- 

CN I line spectrum require a diflFerent explanation. Shock velocities needed to produce ca. 1 

keV -ray in these stars are of order 1000 km/s or so. Such velocity contrasts are diflScult 

to achieve in line-driven instability scenario. 
k>i ■ A more plausible scenario is represented by a Magnetically Confined Wind Shock 

^ ! (MCWS) model, in which stellar wind upflow from opposite footpoints of c losed mag- 



^ _■ netic lo ops collides to form strong, stationary shocks near the loop apex (Babel & Montmerle 



1997bl lah. This model has be e n quit e successful in explaining hard x-rays in Ap and Bp 



stars. However, Gagne et al. (2005) have shown that 2D fully dynamic magnetohydro 
dynamical (MHD) model can also explain both the hardness and the intensity of X-ray 
emission from 6^ Ori C which is an 05.5 star. 

Over the past decade, such 2D MHD models have been applied toward interpreting 
Chandra X-ray spectra of hot stars in general, w hich have resulted in a quite extensive 
series of papers (see ,u d-Doula & Owocki 2003; Iud-Doulall2003l : lOwocki & ud-Doula 



20041 : lud-Doula et al.ll2006l . l2008l .i2009). Although such axisymmetric 2D models work 



well for stars that are slow rotators, they fail to capture more dynamic winds of oblique 
rotators which naturally have - especially when rotation is rapid - lateral structures and 
require full 3D calculations. Here, we present one of the first such fully dynamic 3D 
MHD simulation of a massive star wind, 6^ Ori C. Due to its slow rotation, we are able 
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Figure 1. Snapshot of logarithm of density (left panel) and temperature (right 
panel) taken at time t=400 ks. Note lateral structures are not axisymmetric. Black 
lines represent the magnetic field lines. 



to model its wind as a field-aligned rotator. As we shall see below, even such a slow 
rotator develop lateral structures that otherwise are absent in 2D models. 



2. Chandra X-ray Observations of 6^ Ori C 



Chandra grating observations performed by iGagne et aP (|2005h . have shown that in 
6^ Ori C most of the x-ray emitting plasma is very hot and can be fitted by a global 
temperature of ca. 30 MK. Surprisingly, this temperature does not vary with phase. 
Moreover, the width of spectral lines are very narrow of the order few hundred km/s. 

Bo th of these propert ies are well explained by our fully dynamic 2D MHD simu- 
lations (Gag ne et aDl2005 ). However, the analysis of Chandra data shows that most of 
the hot gas is located around 1.5R:, while our simulations place the hot gas at ~ 2.0R^. 
One of the shortcomings of these 2D simulations is that they impose artificially axisym- 
metry which does not allow for complex structures seen in the simulations to break in 
azimuthal directions which in principle could lead to softer x-ray originating closer to 
the stellar surface. Here we attempt to investigate such a possibility by extending our 
2D simulations to full 3D. 



3. 3D MHD Model 



We study the dynamical competition between field and wind by evolving our 3D MHD 
simulation from an initial condition at time / = 0, when a dipole magnetic field and 
field-aligned rotation is suddenly introduced into a previously relaxed, ID spherically 
symmetric CAK wind. Much of the n u meric al pr ocedures and s tellar parameters used 
are described in ud-Doula & Owockil (l2002h and iGagne et aP (l2005h . However, the 
work presented here in addition includes full energy equations along with radiative 
cooling and a moderate stellar rotation of 10 km/s. We perform our calculation on a 
300x90x90 grid using publicly ava ilable Zeus-MP code, which is the parallel version 
of Zeus-3D dStone & Norman|[l992) . We use spherical geometry. 



3D MHD Simulations 





Figure 2. Slice through the magnetic equator at time t=400 ks showing logarithm 
of density (left panel) and temperature (right panel. 



3.1. Dynamical Evolution 

This artificial, sudden introduction of two dynamically incompatible states - static 
dipole field and outflowing wind - leads in time to formation of certain transient dis- 
continuities in field and flow, e.g., the pointed kinks in field. But these quite quickly 
propagate away within a dynamical timescale of about 25 ks, with the net overall ef- 
fect of stretching the initial dipole field outward, opening up closed magnetic loops 
and eventually forcing the field in the outer wind into a nearly radial orientation. Just 
like in 2D models, within the closed loop regions, wind from opposing hemispheres 
collide along the magnetic equator, shock and eventually cool creating dense equatorial 
structures. Outside and well above the closed region, the flow is quasi steady, although 
now with substantial channeling of material from higher latitudes toward the magnetic 
equator, with max(v6i) > 500 km s"^ even outside the closed loop. This leads to a very 
strong flow compression and thus to a quite narrow equatorial disk of dense, slow 
outflow. But unlike in 2D, after about 100 ks, these structures lose their azimuthal 
symmetry creating instead spike-like structures as seen in Figs. [T]and[2l 

Most of these slow-moving dense structures cannot be driven by radiation. Since 
they are not supported by rotation either, they fall back onto the stellar surface in a 
complex pattern. However, some of the material reach velocities beyond escape speed 
due to momentum deposition in oblique shocks. These dense parcels of material do 
escape the star and may be responsible for X-ray flares observed on some of the massive 
stars. They can also provide a mechanism for creation of clumps, or even be responsible 
for discrete absorption components (DACs). 

3.2. Differential Emission Measure 

In order to understand better the nature of x-ray in our model, we calculate diff'erential 
emission measure (DEM) which essentially is a volume integration of hot gas binned in 
0. 1 dex of temperature. Left panel of Fig. 3 shows the time evolution of DEM for our 
model. Note that there is very little time variation in par with observations. The right 
panel shows averaged over time DEM. Note that most hot gas emit X-ray in the regime 
of 30 MK which again agrees quite well with observations. Further analysis shows that 
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Figure 3. (a) Differential emission measure (DEM) as a function of time. It shows 
that its variation in time is relatively small, (b) Averaged over time DEM for our 3D 
MHD model. Note that the peak DEM temperature ca 30 MK is quite consistent 
with Chandra observations. 



these parcels of very hot gas move very slowly, less than 200 km s"^ suggesting that 
x-ray produced by such gas will form narrow lines. 

Since 3D models allows for mixing in the azimuthal direction, we expected the 
peak temperature of hot gas to be lower than the one obtained from our previous 2D 
MHD models. However, this does not seem to be the case with 3D model suggesting 
even higher temperatures, which agree with the observations even better. Such hot 
gas originate around 0.5-1.0 R:, above the stellar surface, as suggested by the Chandra 
observations. 



4. Conclusion 

Here in this work, we present one of the first fully dynamic 3D MHD simulation of 
the wind of a massive star, 6^ Ori c. Initial analysis suggests that X-ray diagnostics 
from our model in terms of hardness, the location and time variation agree better with 
observations than our previous 2D MHD models. We plan to extend our 3D MHD to 
oblique rotator models, and synthesize dynamical spectrum that can than be directly 
compared with observations. 

Acknowledgments. AuD acknowledges support from NASA grant NNXl 1 AG84G- 
SOl. 
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